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New prediction of extragalactic GeV 7-ray emission from 
radio lobes of young AGN jets 
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ABSTRACT 

We present a new prediction of GeV 7-ray emission from radio lobes of young 
AGN jets. In the previous work of Kino et al. (2007), MeV 7-ray bremsstrahlung 
emission was predicted from young cocoons/radio- lobes in the regime of no coolings. 
In this study, we include cooling effects of bremsstrahlung emission and adiabatic 
loss. With the initial conditions determined by observed young radio lobes, we solve 
a set of equations describing the expanding lobe evolution. Then we find that the 
lobes initially have electron temperature of ~GeV, and they cool down to ~MeV by 
the adiabatic loss. Correspondingly, the lobes initially yield bright bremsstrahlung 
luminosity in ~GeV range and they fade out. We estimate these 7-ray emissions and 
show that nearby young radio lobes could be detected with Fermi Gamma-ray Space 
Telescope. 
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I INTRODUCTION 

Jets in active galactic nuclei (AGNs) are one of the most 
powerful objects in the universe. AGNs with jets are known 
as radio-luminous sources and they have been called as 
radio-loud AGNs. There is increasing evidence that radio- 
loud AGNs significantly interact with ambient medium and 
the cosmological importance of radio-loud AGN feedbacks 
have been advocated (e.g., Best et al. 2007). Among a va- 
riety of radio-loud AGNs, a population with its linear size 
LS < 1 kpc has been called as compact symmetric objects 
(CSOs). The previous studies support a scenario in which 
CSOs propagate from ~ 100 pc scales thrusting away an 
ambient medium and they grow up to Fanaroff-Riley type 

II radio galaxies (FR lis) (e.g., Fanti et al. 1995; Readhead 
et al. 1996a, b; O'Dea & Baum 1997; Stanghellini et al. 1998; 
Snellen et al. 2000; Dallacasa et al. 2000; Orienti et al. 2007). 
CSOs are thus recognized as newly born AGN jets, and they 
are crucial sources to explore physics of radio-loud AGNs in 
their early days. 

On the contrary to various radio observations, high- 
energy emission observations of CSOs had been sporadic 
(Siemiginowska 2008 for review). Recently, XMM-Newton 
observations begin to detect X-ray emissions from some 
of the CSOs although these are not very bright in X-ray 
(Guainazzi et al. 2006; Vink et al. 2006; O'Dea et al. 2006). 
Although a few authors recently indicate possible 7-ray 
emissions from CSOs (e.g., Stawarz et al. 2008), there is little 



observations of CSOs viewed in 7-rays. Hence 7-ray observa- 
tions of CSOs may provide us new independent knowledges 
on the youth era of AGN jets. 

Recently, Kino et al. (2007) (hereafter KKI07) showed 
a possibility of thermal MeV 7-ray bremsstrahlung emis- 
sion from the young radio-loud AGNs. In KKI07, as a 
first step, they focused on the simplest case in which cool- 
ing effects are not significant. These sources roughly cor- 
respond to medium-size symmetric objects (MSOs) with 
1 kpc < LS < 10 kpc and FR II galaxies with LS > 10 kpc. 
Since various cooling timescales in smaller sources tend to 
be shorter than the ones in larger sources, cooling effects 
in CSOs are expected to be more significant than the ones 
in MSOs and FR lis. However little is known about the 
cooling processes in compact radio sources such as CSOs. 
In this study we will consistently solve a set of equations 
which describe young lobe expansions including the effects 
of bremsstrahlung emission and adiabatic loss together with 
initial conditions determined by observations. 

The goals of this paper are (i) to solve a hydrodynam- 
ical evolution of expanding lobes with the cooling effects 
and (ii) to show a new prediction of 7-ray emission from the 
lobes based on the results of (i). In §2, we show a set of basic 
equations which describes an expanding radio lobe inflated 
by its internal energy. In §3, we show the resultant tempera- 
ture evolution of radio lobes. New predictions of bright GeV 
7-ray emissions from young radio lobes is presented in §4. 
Summary and discussion are given in §5. 
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2 MODEL OF RADIO LOBES 
2.1 Basic equations 

Let us consider the dynamics of sideways expansion for 
young radio lobes. We assume that the sideways expansion 
of radio lobe is caused by the internal pressure of the shocked 
jet matter. A breakdown of the internal pressure is an open 
question. The partial pressure of non-thermal electrons alone 
might be small to expand the radio lobes (e.g., Hardcas- 
tle and Worral 2000; Ito et al. 2008 and reference therein). 
Therefore the partial pressure of thermal particles may play 
an important role for lobe expansions. Here we focus on the 
case of thermally expanding lobes and we neglect the partial 
pressure of non-thermal electrons for simplicity. 

The equation of motion and the energy equation, for 
expanding radio lobes, are respectively given by 



(i) 



(2) 



where p cxt R, R — v, P, U t h (= 3P), and Lj are the mass 
density of external ambient medium, the radius, the expan- 
sion velocity, thermal pressure, the thermal energy density of 
relativistic particles of the lobe, and the total kinetic power 
of the jet, respectively. We denote the volume emissivity as 
e coo i and here we focus on the case of e C ooi = et>rcm where 
tbrcm is the bremsstrahlung one. In the case of relativistic 
jets, Lj is given by 



(3) 



where Fj , and Mj are the bulk Lorentz factor of the relativis- 
tic jet, and mass outflow rate of the jet, respectively. The Lj 
is the ultimate source of the lobe expansion. The Lj is a free 
parameter in this model and it is assumed to be constant in 
time. Particular noteworthy compared with KKI07 is that 
we take cooling effects into account. The equation of state 
(EOS) and mass conservation in the lobe are, respectively, 

P = (n- +n+)kT ± + n p kT p , (4) 

^-pR 3 = Mjt, (5) 

where n_, n+, and n p are number densities of thermal elec- 
trons, thermal positrons and total protons in the lobe, re- 
spectively, and p — (n_ + n+)m e + n p m p , T±, and T p are 
total mass density and temperatures of e pairs and pro- 
tons, respectively. Merely for simplicity, we neglect protons 
i.e., n p ~ throughout this work. Although "which plasma 
component is dynamically dominated in AGN jets?" is still 
a matter of debate, previous works at least indicate the ex- 
istence of copious amount of e pairs in jets (Wardle et al. 
1998; Reynolds et al. 1996; Sikora and Madejski 2000; Kino 
and Takahara 2004; Kataoka et al. 2008). Combining the 
above equations all together, we obtain the ordinary differ- 
ential equation of 
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Given the initial values of R, R, R at an initial time to, 
we can solve the evolution of expanding lobes. The initial 
temperature can be obtained via P — p cxt (RR + 3R 2 )/3 
and EOS. 

The free paramter Lj can be normalized as 



L = 5.7 x 10 4 



Mi 



lO- 3 M 



yr- 



erg s 



(8) 



which is typical value for AGN jets. In this work, we will 
examine the range of 5.7 x 10 43 erg s _1 < Lj < 5.7 x 
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10 4D erg s _1 , the corresponding Mj is 1 x 10~ 4 M yr _i < 
Mj < 1 x 10" 1 M Q yr" 1 . The mass outflow rate Mj is the 
essential parameter to determine Lj . 

2.2 Initial conditions of lobes and jets 

Based on recent observations of CSOs (e.g., O'Dea 1998; 
Stanghellini et al 1998; Dallacasa 2000; Snellen et al. 2004), 
we set the initial conditions as follows. Since physical quan- 
tities diverge at R = 0, certain initial conditions of the lobe 
at t — to 7^ are required for solving a set of equations. In 
this work, we treat the evolutions of physical quantities for 
t ^ to and we do not consider t < to- The followings are 
two initial conditions (with subscript of 0) examined in this 
work; 



Ro = 100 pc, 
Ro = 300 pc, 



('() 



Vo 



■■ 0.1c, t = 3.0 x 10 3 yr (Low P ), 
0.2c, t = 4.5 x 10 3 yr (High P ), (9) 



Setting ebrom as for relativistic e pairs (Eq. (22) in Svensson 



and we set dv/dt\t=t ~ 0. This size and velocity is typical 
ones for CSOs (e.g., Stanghellini et al. 1998; Dallacasa et 
al. 2000; Polatidis and Conway 2003; Gugliucci et al. 2005, 
2007; Taylor et al. 2000; Nagai et al. 2006; Kawakatu et al. 
2008). Note that _Ro ~ ^o^o holds here. Since we here deal 
with small radio lobes, the external ambient medium is not 
intracluster medium (ICM) but interstellar medium (ISM). 
By using them, we further set the initial conditions as 

d 2 urr m e P Vo Mjt 

Pa = m p n cxt vo, fei±,o = — : , n_ n = 7^ 77-, (10) 

Mjt 2m e Vo 

where Vo = 4-7r.Ro/3. Hereafter we set the number density 
of external ambient medium n oxt = 1 cm -3 which is typical 
for ISM (e.g., Readhead et al. 1996a and references within). 

In principle, the upper limit of n_ can be constrained 
by the analysis of Faraday depolarization (Dreher et al. 
1987). However, the strong Faraday depolarization observed 
in CSOs (Cotton et al. 2003) are likely to be caused by dense 
foreground matter such as narrow line region. Therefore n_ 
in radio lobes of CSOs has not been clearly constrained. 



3 TEMPERATURE OF RADIO LOBES 

Fig- [U shows the evolution of kT± of expanding lobes which 
is obtained by solving Eq. with the initial condition 
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of the "Low- Pn" case. We examine various outflow rate 



of Mi 



lO" 4 M yr-\ Mi 



10" 



M yr" 1 , Mj 



10~ 2 M yr -1 , and Mj = 10" 1 M yr -1 . Corresponding 
kinetic power of injected jet into the lobes are, respec- 
tively, Lj = 5.7 x 10 42 erg s" 1 , Lj = 5.7 x 10 43 erg 
Lj = 5.7 x 10 44 erg s~~\ and Lj = 5.7 x 10 45 erg s" 1 . The 
initiall temperature of the lobe kT±,o in the case of the Low 
Po case can be analytically evaluated as 



kT±n 



0.18 GeV 



Mj 



lO- 3 M yr- 



PoVo/t 



2.0 x 10 46 erg s~ 



(1 



The resultant kT± evolutions obtained by Eq. © can 
be well approximated by the simple relation of kT± w 
fcTi.o^/to)" 1 + rjm e c 2 /3 for fcT±, > rjm e c 2 /3. In the 
early-phase, the first term (t/to) is dominant for the case 
examined here. Then kT± decreases proportional to 
The agreement of this approximation is due to the negligible 
bremsstrahlung cooling and dominant adiabatic loss in the 
early-phase. We will show the reason in the next section. In 
the late phase, it is found that the temperature asymptoti- 
cally approaches to kT± = Fjm e c 2 /3 = l(iyiO) MeV. This 
clearly coincides with the prediction of "MeV cocoon" by 
KKI07. This implies that the final temperature of lobes is 
determined only by Fj regardless of diverse cooling effects 
in the early phase. Fig. [2] shows the kT± evolution but for 
the "High-Po" case. The initial temperature is given by 



pair expansions. We have neglected the partial pressure of 
proton in Po merely for simplicity. However the neglect does 
not affect the main arguments in this work. Because the pro- 
tons also simply cool down by the adiabatic loss for t > to- 
To form "High-Po" lobe, very fast Tj ~ 50 — 100 is required 
when assuming that hot protons with kT v = FjiripC 2 /3 are 
the heating source of e pairs. The required Fj is by a fac- 
tor of ~ 2 — 3 larger than the fastest Tj ~ 30 estimated 
from VLBI observations (Kellermann et al. 2004) and it is 
comparable to Tj ~ 50 — 100 indicated by rapid variabilities 
of blazars (e.g., Begelman et al. 2008; Ghisellini and Tavec- 
l)chio 2008). Since we take the oversimplified assumption of 
proton neglect, it is not possible any further to explore the 
hot lobe formation at t < to. We keep this as a subject for 
future investigation. It is worth to mention that youngest 
radio sources termed "high frequency peakers" (HFPs) with 
ages of < 10 3 yr have been recently observed (e.g., Orienti 
et al. 2007) and HFPs may give us substantial hints for the 
physics in t ^ to- 



4 7-RAY EMISSION FROM RADIO LOBES 

Using Eqs. (|10|l . (|11[) , and (|12[) , we can estimate the bolomet- 
ric bremsstrahlung luminosity Lbrcm.o = £brem,oVb at t = to 
for the "Low Po" case 



kT ± , = 13 GeV 



M s 



lO" 3 M yr- 



PoVo/to 



1.4 x 10 48 erg s" 1 



(12) 



The behavior is essentially the same as Fig.[T] Adiabatic loss 
phases for the "High-Po" case last longer than those for the 
"Low-Po" case simply because of its higher kT± t o. 

Next, we show how kT± t o will change when we choose 
different initial conditions. The terms Mj and PoVo/to are 
the only ingredients which determine the initial tempera- 
ture kT±.o- The predicted feT±,o will increase by the increase 
of Po, «o, and p cxt according to PoVo/to oc p cxt RoVo, and 
vice versa. Therefore, in principle, kT± t o can become higher 
or lower than MeV/GeV ranges examined in Fig. [T] When 
kT± t o is order of keV, the predicted kT± t o is smaller than 
rjm e c 2 /3 = l(iyi0) MeV. Therefore, the lobe plasma is 
gradually heated by the jet injection and kT± asymptoti- 
cally approaches to Fjm e c 2 /3 = l(Fj/10) MeV. However, 
the case of Mj ~ 10 M yr -1 requires Lj ~ 10 48 erg s _1 
(see Eq. (|SJ) which tend to exceed the ones for very power- 
ful AGN jets (e.g., Rawlings and Saunders 1991; Ito et al. 
2008). Here we do not discuss this case further. 

Lastly we comment on the formation of initial lobes. In 
order to accomplish the initial hot lobes at t = to, a heat- 
ing process for e pairs is needed at t ^ to. Theoretically, 
a mechanism of electron heating (and acceleration) in a col- 
lisionless shock is a matter of debate (e.g., Shimada and 
Hoshino 2000; Ohira and Takahara 2007). If pairs are effec- 
tively heated up by various kind of plasma instabilities via 
hotter protons (e.g., Lyubarsky 2006), it could be possible 
to attain kT±.o ~ kT Pt o. It is beyond the scope of this paper 
to explore microscopic processes in the plasma. Instead, we 
have uniquely determined T±,o by using the quantities re- 
quired from the observations and the assumption of thermal 



ibrom.O = 1-5 X 10 
kT±n 



(of^v)[' + '»( 



2.7 x lO" 2 cm- 3 
1 



Hi 



Vo 



2 x 10 B2 cm 3 , 
erg s \13) 



with the same M — 10 3 M yr 1 have hotter e ± pairs and 
smaller n_,o in larger Vo. Then, the "High Po"case leads to 
brighter L brO m,0 of 



-^brcm.O — 



1.4 x 10 

kT ± , 



(i 



5 x lO" 3 
l+ln| 0.10 



V 



3.3 x 10 63 



kT- 



±,o 



13 GeV 



erg s 



(14) 



Fig. E3 shows the bremsstrahlung spectra at t = to. Adopted 
Mj are same as in Fig. 1. The source is located at the 
distance of 10 2 Mpc which corresponds to nearby ob- 
served CSO samples (e.g., Snellen et al. 2004). For these 
cases, the predicted spectra are brighter than the sensitivity 
of Fermi/LAT (http://www-glast.stanford.edu/). In other 
words, we predict radio lobes in CSOs as a new popula- 
tion of GeV-7 emitter in the universe. From Eqs. (113|) and 
(|14[) . we find that the increase of Mj leads to the enhance- 
ment of I/brcm,o since I/brem.o oc Mj holds where we neglect 
the logalithm terms in (|13[) and (|14fl . We can also say that 
lower kT±fi leads to brighter Lbrcm,o because n± i0 oc Mf 
dominates kT±,o oc Mr 1 , and vice versa. Regarding the 
luminosity evolution, we obtain Lbrcm oc kT±n±V oc n± 
where we use the relation of kT± oc t~ and nV oc t for 
the early phase. We thus find that Lbrem decrease with time 
and it will fade out because the lobe becomes dilute. This 
sort of negative luminosity evolution has been well known 
for synchrotron emission (Readhead et al. 1996b; Begelman 
1996) and bremsstrahlung emission (KKI07). Once we ob- 
tain Lbrom.o, we can easily check that the radiative cooling 
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Figure 1. Evolution of feT± for the "Low-Po" case which 
starts from the age of to = 3000 yr. The examined cases are 
A/j/(Mq yr" 1 ) = 10" 4 , IO" 3 , IO" 2 , and 10" 1 and they are 
labeled as L4, L3, L2, and LI, respectively. Corresponding ki- 
netic power of the jet are, respectively, Lj = 5.7 X 10 42 erg s _1 , 
Lj = 5.7 x 10 43 erg s" 1 , Lj = 5.7 X 10 44 erg s" 1 , and Lj = 
5.7 X 10 45 erg s — 1 . The temperature decrease by the adiabatic 
loss are seen in the initial phase. Each case asymptotically goes to 
constant kT± phase which is predicted by KKI07. Bremsstrahlung 
cooling is found to be negligible. 
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Figure 2. Evolution of kT± for the "High-Po" case which 
starts from the age of to = 4500 yr. We examine the cases 
of Afj/(M yr" 1 ) = 10" 4 , IO" 3 , IO" 2 , and IO" 1 and they 
are labeled as H4, H3, H2, and HI, respectively. Correspond- 
ingly, kinetic power of injected jet into the lobes are, respec- 
tively, Lj = 5.7 X 10 42 ergs -1 , Lj = 5.7 X 10 43 ergs -1 , 
Lj = 5.7 X 10 44 erg s" 1 , and Lj = 5.7 X 10 45 erg s" 1 . Similar 
to the "Low-Po" case, adiabatic loss is dominant in the initial 
phase and bremsstrahlung cooling is negligible. 



is effective or not. Since Lbrem decreases with time, the ef- 
fect of bremsstrahlung is most significant at t — to- At the 
time, Lbrem.o *C Lj holds where Lj is the source term of 
the internal energy flux injected in the lobe. Therefore we 
find that back-reaction of bremsstrahlung cooling is negligi- 
ble and the behavior of kT± in Figs. [1] and [2] are governed 
by the adiabatic loss. 

Lastly we discuss the change of Lb re m,o for different 
choices of initial conditions. Two ingredients to determine 
ibrom.o are n-fi and kT± t o for given Vo. Concerning kT±,o oc 
PoVo/to, the predicted L hrcul increases as Ro, vo, and p ex t 
becomes larger according to kT±,o oc pcxt-Ro^O) an d vice 
versa. Hence it is worth to note that lobes with faster side- 
ways expansions will shine brighter. Indeed there are some 
CSOs with Vo > 0.1c (Polatidis and Conway 2003; Gugli- 
ucci et al. 2005, 2007). The quantity n_,o is the important 
parameter, since n_,o dependence of Lbrcm.o is significant as 
ibrom.o oc n 2 _o- For the estimate of ri-,o, we assumed the e 
pair jet so far. If instead we assume the pure electron/proton 
jet for given Lj and Fj, then n_,o decreases by a factor of 
m e /m v . Then L br em,o decreases by a factor of (m e /m p ) 2 , 
and the predicted Lbrcm.o becomes much smaller than the 
Fermi/LAT sensitivity. 



5 SUMMARY AND DISCUSSION 

We have investigated the temperature and luminosity evo- 
lutions of radio lobes of CSOs which expands by their own 
thermal pressure. In the previous work of KKI07, MeV 7-ray 
bremsstrahlung emission was predicted from young lobes in 
the regime of no coolings. In this work, we include cooling 
effects of bremsstrahlung emission and adiabatic loss. Below 
we summarize the main results of the present work. 

(i) We examine the evolution of kT± together with the 
initial conditions determined by observed CSOs. By solving 




hv [eV] 



Figure 3. Predicted bremsstrahlung emission from young radio 
lobes at t = to with the distance 10 2 Mpc. The spectra for "High- 
Pq " case (HI, H2, H3, and H4) are brighter than the sensitivity 
of Fermi Gamma-ray Space Telescope while the spectra for "Low- 
Po" case (LI, L2, L3, and L4) are less luminous for the detection. 

a set of equations describing the expanding lobe evolution, 
it is found that the lobes initially have electron temperature 
of ~GeV for Fj ~ 10, the lobes then cool down to MeV by 
the adiabatic loss. During the early phase, kT± is governed 
by the adiabatic loss alone. Since the adiabatic loss is more 
effective than bremsstrahlung cooling in any case. In the 
late phase, kT± asymptotically approaches to the constant 
temperature of Tjm e c 2 /3 = 1 (Tj/10) MeV which has been 
predicted in KKI07. 

(ii) Thermal bremsstrahlung emission peaked about 
GeV-7 band is newly predicted in CSOs because n± and 
kT± in younger radio lobes are larger than those in older 
ones. These spectra can be brighter than the sensitivity of 
Fermi/LAT for nearby CSOs in the case of High-Po lobe. 
This means that young radio-lobes can be a new population 
of GeV-7 emitter in the universe. From Eq. (|13[) showing 
Lbrom.o, we see that the radio lobes with larger vo (i.e., faster 
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expansions) and/or the one with larger n_ (i.e., larger Lj) 
yield brighter emission for given Ro- As for a specific source, 
we stress the importance of young recurrent lobe 3C 84 (e.g., 
Asada et al. 2006). Since 3C 84 is located at the center of 
very nearby Perceus cluster (z = 0.018), a deep observation 
of the Fermi/LAT collabolated with radio obserbations on 
it will give us tight constraints on the physics of young radio 
lobes. Near future mission of the VLBI Space Observatory 
Programme 2 (VSOP-2) with unprecedented high angular 
resolution (Tsuboi et al. 2009) would play important roles 
for direct constraint on vo and Ro of HFPs. 

(iii) Stawarz et al. (2008) investigated a variety of non- 
thermal 7-ray emissions from lobes of CSOs. They claim that 
the predicted non-thermal emissions can be also detected 
by Fermi/LAT with an assumed certain electron accelera- 
tion efficiency. In GeV 7 range, the dominant components 
are inverse-Compton scattered emissions of ultraviolet and 
infrared photons. It is worth to note that even though the 
luminosity of thermal bremsstrahlung and inverse-Compton 
ones are comparable, the spectrum shape thermal compo- 
nent is quite different from non-thermal spectra. Hence it is 
straightforward to distinguish whether the emission is ther- 
mal or non-thermal one. 

(iv) The observations of Fermi/LAT will be tests for some 
unresolved questions of AGN jets. Suppose the case that we 
exactly know Ro, vo, and n oxt , and a source distance. If we 
observe the predicted thermal GeV-7 emission from CSOs, 
then it suppose the scenario in which CSOs have relativistic 
jets and their lobes thermally expand. If we do not detect 
it, it is attribute to lower kT± and/or smaller n±. Possible 
reasons are as follows; (1) the jet is mainly made of e/p 
plasma with the same Lj, (2) the jet consists of e ± plasma 
on the whole but with smaller Lj , (3) the lobe is expanded by 
energetic non-thermal particle and the actual n± is smaller, 
(4) the jet has non-relativistic speed which leads to the lower 
kT±, (5) other cooling processes could make kT± lower. We 
remain them as our future investigations. 
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